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Context. Photometry and polarimetry have been extensively used as a diagnostic tool for characterizing the activity of comets when 
they approach the Sun, the surface structure of asteroids, Kuiper-Belt objects, and, more rarely, cometary nuclei. 
Aims. 133P/Elst-Pizarro is an object that has been described as either an active asteroid or a cometary object in the main asteroid 
belt. Here we present a photometric and polarimetric study of this object in an attempt to infer additional information about its origin. 
Methods. With the FORS1 instrument of the ESO VLT, we have performed during the 2007 apparition of 133P/Elst-Pizarro quasi- 
simultaneous photometry and polarimetry of its nucleus at nine epochs in the phase angle range ~ 0° - 20°. For each observing epoch, 
we also combined all available frames to obtain a deep image of the object, to seek signatures of weak cometary activity. Polarimetric 
data were analysed by means of a novel physical interference modelling. 

Results. The object brightness was found to be highly variable over timescales < 1 h, a result fully consistent with previous studies. 
Using the albedo-polarization relationships for asteroids and our photometric results, we found for our target an albedo of about 
' 0.06-0.07 and a mean radius of about 1.6 km. Throughout the observing epochs, our deep imaging of the comet detects a tail and an 

anti-tail. Their temporal variations are consistent with an activity profile starting around mid May 2007 of minimum duration of four 
months. Our images show marginal evidence of a coma around the nucleus. The overall light scattering behaviour (photometry and 
polarimetry) resembles most closely that of F-type asteroids. 

■ Key words. Comets: individual: Comet 133P/Elst-Pizarro - Polarization - Scattering 
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' 1. Introduction trinsic brightness may increase above the normal linear brighten- 
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ing with phase angle. Both phenomena are a consequence of the 



^ . Light scattering of minor solar system bodies, such as asteroids, surface micro structure and usually attri buted to the combined 

O comets, and Kuiper-Belt objects, plays a central role in deter- action of the shadowing e ffect and the coherent backscattering 

O : mining global physical parameters such as size and albedo, and of light At larg£r phase angle valueSj the yalue of me polar . 

. the detailed understanding of the surface structure such as mi- ization minimumi the inversion angle at which the polarization 

L; ; cro structure and single scattering albedo of the body surface changes f rom being parallel to the scattering plane to becoming 

. , _h . d Muinonen || 2004|) . Thls a PP lles t0 both remote sensin § °bserva- per pendicular to it, and the slope of the polarimetric curve at the 

^ ■ tions from Earth and to in situ measurements during spacecraft cross . over point differ f or asteroids of different surface taxon- 

H ; encounters. Light scattering is studied by means of photomet- omy and may even be used to assess the albedo of the object 
»t i ric and linear polarimetric measurements obtained at different 

phase angles, i.e., the angle between the Sun, the object, and the The systematic research and classification of the linear po- 
observer. The way in which photometry and linear polarization larimetric and photometric phase functions for the zoo of small 
change as a function of the phase angle helps us to characterize bodies in the solar system remains in its infancy, mostly be- 
the scattering medium, in both the case of solid surfaces and dust caU se of a lack of good coverage by measurements for the 
ejecta around comets. various objects types to be considered (for instance, asteroids, 
The phase angle range between 0° and about 30°is of particu- cometary nuclei, Kuiper-Belt objects, Trojans, minor satellites), 
lar interest. At small phase angles (^ 10°), the linear polarization Photometric and pol arimetric techniques h ave been used to 
of small bodies in the solar system is usually parallel to the scat- study aste roids (e.g.. iMuinonen et alJ 120021). Kuiper-Belt ob - 
tering plane, in contrast to what is expected from the simple sin- jects (e.g., iBoehnhardt et al.l l2004; Bagnulo et ai1 l2006l [2008), 
gle Rayleigh-scattering or Fresnel-reflection model, and the in- and the activity of comets while approaching the Sun, (e.g. 
iPenttila et al.l 120051) . In contrast, photometric and polarimetric 



* Based on observations made with ESO Telescopes at the Paranal data of cometary nuclei are scarce in the literature. Photometric 
Observatory under programme ID 079.C-0653 (PI=Tozzi) phase functions of cometary nuclei are available for fewer than a 
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dozen of objects (see Lamy et alJ 2004: 


Boehnhardt et al. 


2008; 


Lamv et al.ll2007aLbUch iLi et al.ll2007alb 


; Snodgrass et al. 


2008; 


Tubiana et al.l2008). The polarization of cometary nuclei is even 



less observed and studied than photometry. To the best of our 
knowledge, only a single comet nucleus, that of 2P/Enck e, has 
been observed in polarized light (Boehnha rdt et al.l l2b08). over 
an extented phase-angle range. Although the shape of the po- 
larimetric curve of 2P/Encke resembles to those of asteroids, the 
numerical values of the polarization minimum, slope, and in- 
version phase-angle differ remarkably from those of other small 
bodies (including asteroids) in the solar system measured so far. 

Here we present new observations of 133P/Elst-Pizarro, an 
object classified as both a comet and an asteroid (7968 Elst- 
Pizarro) and seems to be by its nature i n between comets 
and asteroids. It is a small km-size body dHsieh et alJ [2009) 
that belongs to the ty pe of so-called main-belt comets (MBC, 
iHsieh & Jewittl 120061) of which only five objects are known so 
far - P/2005 Ul (Read), 176P/( 118401) LINEAR, P/2008 Rl 
(Garradd), and P/2010 Al (LINEAR) are the four others. It has 
an orbit rather typical of main-belt asteroids, for which reason 
it is also designated as asteroid, is similar to the Themis col- 
lision family, or a s ub-family of it, called the Beag le family, 
dNesvornv et all 2008). It was discovered in 1996 by lElst et alJ 
(1996) at the European Southern Observatory (ESO) in La Silla 
(Chile); at that time it displayed a long thin dust tail, which trig- 
gered its initial designation as a comet. Although details of the 
activity's origin are as yet unknown, the activity producing the 
tail seems to be recurrent and occur close to perihelion passage, 
while the object m ay be inactive over the remainder o f its or- 
bit around the Sun dHsieh et alJl2004tlJewitt et alj|2007l) . Rather 
than a mostly inactive comet, 133P/Elst-Pizarro is nowadays 
more commonly regarded as an asteroid with an ice reservoir, 
whic h periodically sublimates, with consequent material ejec- 
tion dBoehnhardt et al.lll997l:lHsieh et alJl2004tlToth Il2006h . 

2. Observations 

Polarimetric and photometric observations of 133P/Elst-Pizarro 
were obtained in service mode at the ESO Very Large Telescope 
(VLT ) from May to Septemb er 2007, with the FORS1 instru- 
ment ( AppenzeUeLeTalll 998), using broadband Bessell R and V 
filters. Linear polarization measurements were obtained at nine 
different epochs, one before perihelion (in May 2007), and the 
remaining eight after perihelion (from July to September 2007). 
Each series of polarimetric observations consisted of a 30 s ac- 
quisition image obtained in the R Bessell filter without polari- 
metric optics, and a series of images obtained with the half wave- 
plate set at 12 to 24 position angles in the range 0°- 337.5°, 
in steps of 22.5°, both with R and V Bessell filters. We also 
obtained photometric imaging in R and V Bessell filters (i.e., 
without polarimetric optics) at the same nine epochs (quasi- 
simultaneous to the polarimetric series). 

We performed a preliminar inspection of the Line of Sight 
Sky Absorption MonitoiQ (LOSSAM) plots of each observing 
night showing the atmospheric conditions on site, and we found 
that for all epochs, except one, sky transparency was close to 
photometric at the time of the observations. Night 29 to 30 
August 2007, corresponding to our target at phase-angle = 15.3°, 
was cloudy. For polarimetry, the only impact is in terms of a 
reduced signal-to-noise ratio, while photometric measurements 
and tail length measurements obtained during that night should 
not be considered reliable. Some of the frames obtained on 24 



September could not be used because of background objects 
overlapping the image of 133P/Elst-Pizarro. 

3. Data analysis 

Both polarimetric and photometric data, including acquisition 
frames and science frames, were pre-processed in a similar way. 
Frames were bias-subtracted using a master bias obtained from 
a series of five frames taken the morning after the observations, 
then divided by a flat-field obtained by combining four sky flat 
images taken during twilight with no use of polarimetric optics. 
From this point, we proceeded with different strategies tailored 
to the specific analysis that we aimed to perform. 

3. 1 . Deep imaging and tail analysis 

Deep imaging of 133P/Elst-Pizarro was obtained by coadding 
the polarimetric images as follows. For each frame, we consid- 
ered separately the two images with opposite polarization that 
are split by the Wollaston prism, from which we subtracted the 
background using SEXTRACTOF0 with the mesh option. "Full res- 
olution" background images were checked to ensure that they 
were constant around the target, and that background subtraction 
neither remove nor added small scale features. We then coad- 
ded the two images with opposite polarization for each of the 
frames obtained at various positions of the retarder waveplate. 
These frames were corrected for their corresponding airmass 
values. Extinction coefficients, colour correction terms, and ze- 
ropoints were obtained from the observations performed within 
the framework of the FORS1 instrument calibration plan. For 
the R filter, we adopted the values of k R = 0.087 + 0.014 and 
ky-R = -0.097 + 0.01 for the extinction and the colour correc- 
tion, respectively; for the V filter, we adopted ky — 0.1 12+0.012 
and k B -v = -0.015 + 0.007. For the colour indices of 133P/Elst- 
Pizarro we adopted V - R — 0.41 and B — V — 0.60, which are 
close to solar colours. However, the instrument calibration plan 
does not include photometric calibration of the images obtained 
with the polarimetric optics, and the transmission function of 
the polarimetric optics is documented in neither the FORS user 
manual nor the literature. We therefore performed a quick pho- 
tometric calibration of the images obtained with the polarimetric 
optics by comparing the flux of the same objects observed with 
and without polarimetric optics. We estimated that the absorp- 
tion of the polarimetric optics is about 0.60 and 0.54 mag in R, 
and V, respectively. Instrument zeropoints for the polarimetric 
mode were thus calculated by subtracting 0.60 or 0.54 from the 
night zeropoints obtained within the context of the instrument 
calibration plan in imaging mode. We note that this correction 
was obtained without accounting for a change in the extinction 
and colour coefficients caused by the polarimetric optics. The 
impact of this approximation should be minimal, relative to the 
uncertainties discussed in Sect. 13.31 

All frames obtained during a single epoch were average- 
combined using a cr-clipping algorithm, adopting the pixel me- 
dian value as a center for clipping, and calculating the average of 
the non-clipped pixels. These combined images were then cali- 
brated in Af, where A is the albedo and / the filling factor, using 
Ea. mof lTozzi etal.ld2007h 

Af - 5.34 x 10 11 (r/d p f 10 x (m ^ zp "' ) N , (1) 

where r is the heliocentric distance in AU, m is the absolute 
magnitude of the Sun in the considered filter, ZP,„ is the zero- 



1 Online data available at the ESO web site 



2 http://astromatic.iap.fr/software/sextractor 
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Fig. 2. Contour plots of Af in Bessell filter for the combined 
image obtained on July 17. Contour levels correspond to Af = 
2 x 1(T 9 , 4 x 1(T 9 , and 6 x 1(T 9 . 
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Fig. 1. Contour plots of A/ in /? and V Bessell filter (left and 
right columns, respectively). North is up, and east to the left. 
The black arrow and red arrow represent the direction of the Sun, 
and the negative target velocity as seen from the observer in the 
plane of the sky, respectively. The x and y axes are in units of 
10 3 km. Contour levels are calibrated in Af (for explanations see 
text). Contour levels correspond to Af = 7.0 x 1(T 9 , 8.3 x 10" 9 , 
1.0 x 10~ 8 , 1.25 x 10~ 8 , 1.65 x 10~ 8 , 2.5 x 10~ 8 , and 5 x 10~ 8 . 



point in that filter, and d p is the CCD pixel scale in arcsec (0.25" 
in our case), and N is the number of electrons per pixel. 
Images were finally magnified using a scale factor 



s = 0.725 A d D 



(2) 



where A is the geocentric distance in AU at the time of the spe- 
cific observation. In this way, the pixel coordinates where con- 
verted into identical spatial coordinates expressed in 1000 km. 

FigureQ]shows the contour plots of the coadded polarimetric 
images of 133P/Elst-Pizarro at the nine observing epochs, i.e., 
from 38 days before to 87 days after perihelion. 

The image of 22 May 2007 shows a marginal indication of a 
tail in the anti-solar direction. In all images obtained after perihe- 
lion, we clearly detected one or two narrow tails, either at posi- 
tion angles PA corresponding to the direction of the Sun, at 180° 
from it, or in both directions. Since the tails are seen in both V 
and R filters, we conclude that they represent dust-reflected sun- 
light (and not so much - or possibly not at all - from gas emis- 
sion). Table[T]provides geometric information about the two tail 
features measured from the images. 

For comparison purpose, in Fig.[T]we adopted the same con- 
tour levels for all epochs, 7 x 10~ 9 being the one at the small- 
est Af value. However, we note that in all images, apart from 
one obtained before perihelion, tails extend to a greater distance 
than indicated by the 7 x 10~ 9 contour level (although, below this 
value, background noise becomes quite significant). The tail was 
detected with the highest signal-to-noise ratio in the two images 
obtained in July 2007, extending at least up to 25 000 km from 
the photometric centre of the object, pointing toward the Sun. 
We note that because of projection foreshortening, the measured 
length of the tail in the sky compares to a much longer extension 
in space. In the images obtained on July 13, and in the regions 
at projected distance between 6 000 and 20 000 km, we estimate 
that the scattering cross-section of the tail (i.e., the projected sur- 
face multiplied by the albedo) is about 0.27 km 2 and 0. 19 km 2 in 
the R and V filter, respectively. In the images obtained on July 17 
for the same regions, we measure about 0.24 km 2 and 0.21 km 2 
in the R and V filter, respectively. Figure [2] shows the contour 
plot for the combined image obtained on July 17 in the R filter 
setting for the lower contour level curve of value Af = 2 x 10~ 9 . 

In August, we detected, apart from the primary tail, a weak 
secondary tail, which then prevailed in brightness over the first 
tail in the two images obtained in September 2007. Hereafter, 
the primary one will be referred to as "Tail 1", the secondary 
tail as "Tail 2". Tail 1 points westward, i.e., close to opposi- 
tion it is directed towards the Sun. It appears as an anti-tail (a 
sun-ward pointing dust tail) in our July 2007 image, and there- 
after as a normal tail. It is brighter than Tail 2 before about 
mid-August 2008, then fades away and is no longer detectable 
in our last exposure series on 24 September 2007. Tail 2 in the 
eastern hemisphere always appears as an anti-tail. It is first im- 
aged in early August 2007 as a short eastward extension in the 
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Table 1. Geometric information about the dust tails of 133P/Elst-Pizarro. 





Phase 






Tail 1 

ran 1 




Tail 1 

ran z 


Date 


Angle 


OA (1) 


PA' ' 


* j l .0,(3) 

projected length' ' 


PA' ' 


projected length' ' 


yyyy mm dd 


C) 


(°) 


(°) 


(10 (3) km) 


o 


(10 (3) km) 


2007 05 22 


19.7 


256.2 


257 


- 


- 


- 


2007 07 13 


3.4 


248.1 


255 


18 






2007 07 17 


1.7 


237.7 


256 


17 






2007 08 01 


5.0 


85.7 


257 


12 






2007 08 05 


6.7 


84.2 


260 


12 


81 


3 


2007 08 09 


8.2 


83.4 


258 


12 


82 


4 


2007 08 30 


15.3 


81.9 


262 


7 


82 


6 


2007 09 19 


19.6 


81.0 


261 


4 


81 


10 


2007 09 24 


20.3 


80.7 


261 




81 


10 



(1) Position angle of the antisolar vector. (2) Position angle of the tail. (3) Approximate tail extension obtained using the contour level 
Af = 7 10~ 9 of Fig.[T| averaged in the R and V filters. 



isophote pattern in Fig. \T\ then brightens above Tail 1 by the 
end of August, and remains detectable until the end of our ob- 
servations in late September 2007. The appearance of the two 
tails in 133P/Elst-Pizarro resembles in terms of its behaviour the 
du st phenomena observed in this object in 1995 and described 

bv lBoehnhardt et all d!997l) . 

U sing the Finso n-Probstein (FP) code dBeisser & Drechsefl 

Il992t lBeisser|[l990l) . we could develop a qualitative and in part 
also quantitative understanding of the dust activity of 133P/Elst- 
Pizarro. The FP code allows one to calculate and display the so- 
called synchrone-syndyne pattern of cometary dust tails, where 
synchrones represent the location of the dust emitted from the 
nucleus at the same time and syndynes represent the dust sub- 
ject to the same solar radiation pressure characterized by the fi r 
parameter, the ratio of the force of solar radiation pressure to that 
of the gravity of the Sun. During perihelion passage in 2007, the 
synchrone-syndyne pattern of 133P/Elst-Pizarro shows a very 
dense and narrow grid of overlapping lines, indicating a very 
narrow and spiky dust tail. It also suggests that the two observed 
tails are phenomena of projection caused by a single dust tail, 
formed by the dust activity of the object that lasted over a cer- 
tain period around perihelion. Owing to the low orbital inclina- 
tion of 133P/Elst-Pizarro, the dust tail is seen almost edge-on 
from Earth during the whole observing interval. Depending on 
the observing dates, the FP calculations show that the two tails 
represent different dust grain populations, emitted by the nucleus 
at different time before the observing epoch and located on dif- 
ferent sides of the nucleus as projected onto the sky of the ob- 
server on Earth. The projection-induced transition between Tails 
1 and 2 happened shortly after opposition passage of the object 
on 20 July 2007. 

From the images, we conclude that the dust activity of 
133P/Elst-Pizarro was very low (or even absent) during our first 
observing night (22 May 2007) 38 days before perihelion (top 
panels of Fig.[TJ. However, the presence of Tail 1 on 30 August 
2007, and possibly on 19 Sept. 2007, suggests that significant 
dust activity had started shortly (a few days) after 22 May 2007, 
since the dust grains seen in this tail region must have been re- 
leased by the nucleus about 100 days or more before the date of 
observation, i.e., in late May 2007. On the other side, the exis- 
tence of Tail 2 in September 2007 is indicative of dust release 
by the nucleus that was still ongoing by the time of our observa- 
tions, i.e., almost 90 days after perihelion passage. 

The extension of Tail 2 from the nucleus allows us to esti- 
mate an approximate maximum y6 r value for grains in the most 
distant region of this tail: we found a maximum /3 r of about 0. 15, 



which is characteristic for instance of a few fim or ten fim size 
grains from silicate or absorptive materials. On 9 August 2007, 
the width of Tail 1 is about 1000 km at a projected distance of 
5000 km from the nucleus indicating a slow (out-of-plane) ex- 
pansion speed of the dust of only about 1 .5 m/s. The tail width on 
13 and 17 July 2007 at the same distance range (4000-5000 km) 
was measured to be 1750 and 2250 km, respectively, which re- 
sults in an expansion speed out of the orbital plane of about 1 .45 
ms _1 . Both results support the conclusions of a low expansion 
velocity of the dust grains in this object. 

Finally, to obtain tighter constraints of the tail(s) brightness, 
we calculated the average Af at different distances along the di- 
rection identified by the tail(s). We measured the average flux in 
rectangular areas of 20 pixels, five along the direction identified 
by the tail, and four in the direction perpendicular to it0 Figure|3] 
shows the Af in the R and V filters, and the colour index V — R, 
versus distance from the photometric centre of the object. The 
points at positive distances refer to Tail 1, and those at negative 
distances refer to Tail 2. The V — R points of the rightmost pan- 
els were plotted only when flux was detected in both filters at the 
minimum 2 cr level. The V — R colour of the tails is about 0.4. 
It is mostly constant along the tail axes, and close to a neutral 
intrinsic colour typical of a flat spectrum in the visible. It is also 
in good agreement with the V — R colour of the nucleus itself 
(see Sect. 13.3b . The neutral colour of the dust tails is compat- 
ible with light scattering by grains that were much larger than 
the wavelengths of the filter measurements. It also implies that 
the surface material and the dust released by 133P/Elst-Pizarro 
is not "red" as frequently observed for cometary dust and con- 
sidered indicative of the space weathering effects on the surface 
materials. 

3.2. Searching for coma activity 

We provide a formal description of our approach to evaluate 
coma contribution from our mea surements. This is based on a 
strategy originally developed by Toz zi & Licandrol d2002l) and 
Tozz i et alj d2004l) . Thereafter, we describe the results of the 
analysis performed on the combined polarimetric images of 
133P/Elst-Pizan-o. 

We assume that the number of detected electrons e~ of the 
object per unit of time within a circular aperture of radius a, 
Niot(a), is the sum of the contribution of the nucleus, A/n, plus 



3 We note that the measured tail Af values may be possibly overesti- 
mated by up to 10 % because of a faint coma (see Sect. 13. 2\ . 
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Fig. 3. Tail Af measured along the direction identified by the tail, with respect to the photometric center, and tail colour indices 
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panels, they correspond to the colour index derived for the nucleus. 



the contribution of the coma, Nc and, possibly, a spurious con- 
tribution A/b, due to non-perfect background subtraction 

NUa) = N n (a) + N c {a) + N z (a) , (3) 

where a is the rad ius of the aperture in pixels. Following 
lA'Hearn et a fl (fl98l . the flux of a (weak) coma around the nu- 
cleus in a certain wavelength band can be written as 

Fc =M^)>- (4) 

where A is the mean albedo (unitless), / is the filling factor (unit- 
less), A is the geocentric distance and p is the projected distance 
from the nucleus (corresponding to the aperture in which the 
flux Fq was measured), expressed in identical units, r is the he- 
liocentric distance expressed in AU, F Q is the solar flux at 1 AU 



integrated in the same band as Fq, and convolved with the filter 
transmission curve. Equation (0 applies to the reflected light of 
the dust coma, which in the visible clearly dominates over contri- 
butions from gas emission. Based on the hypoth esis of uniform 
and is otropic ejection of dust at constant velocity. [A'Hearn et al.l 
(Il984l) note that the product Afp is constant with p. The back- 
ground is assumed to be constant in the small region measured 
around the object, hence, linearly proportional to the aperture 
area. Equation (01 can thus be written 

ATwCa) = N N (a) + k (C) a + k (B) a 2 , (5) 

where k (C) and £ (B) are terms independent of a. We then trivially 
obtain 

KM) = iKoM) = %{a) + + 2ak™ , 

(6) 

KM) = faNnffl) = N"(a) + 2kW . 
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Plotting N^ t (a) versus a allows one to estimate whether back- 
ground subtraction is optimum, as at large apertures compared 
to the comet image - barring the presence of background ob- 
jects - A/^ t (a) should converge toward the £ (B) value, which for 
optimum sky subtraction should be zero. 

Assuming that the nucleus is a point source, the term N^(a) 
is the point spread function, and, for apertures sizes that are large 
compared to the seeing, should tend to zero. If an extended coma 
is present, we expect a contribution from the term N' c {a) con- 
stant with a. The coma contribution can then be evaluated by 
measuring the constant & (C) as a weighted average of N' lol {a) in 
the aperture interval [a,,, a,- J, with a,-, < a,,. The k (C) value can 
finally be used in Eq. © to distinguish, for each aperture value, 
the flux due to the coma and that due to the nucleus. In particu- 
lar, the former one increases linearly with aperture, and the latter 
should appear constant at the aperture values at which N{ ot {a) is 
constant, and which were used to determine k (C K 

The nucleus contribution An can then be transformed to 
magnitude m, following the standard recipes for photometry. 
The coma contribution Ac can be transformed into astrophys- 
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Fig. 5. Afp for the coma on 22 May 2007. Red circles refer to R 
filter, blue squares to V filter. Dashed lines mark the average Afp 
values in the two filters, obtained neglecting the points marked 
with empty symbols. 



ically meaningful te rms using the quantity Afp introduced by 
lA'Hearn et ail(ll984h . using the formula 



Afp = 1.234 x 10 19 10 



f>A(m -TP m ) r l 



4 I £ (C) ■ 



(7) 



where r and A are measured in AU, d p in arcsec per pixel, £ (C) in 
e~ per pixel, and Afp is obtained in cm. 

To test this algorithm, we used the frames obtained on 22 
May 2007, when 133P/Elst-Pizarro exhibited the least evidence 
of tail activity, and its image appeared relatively isolated from 
background objects. We applied our algorithm to both 1 33P/Elst- 
Pizarro and an isolated background star of similar brightness. 
Suitable comparison objects are extremely rare in our images. 
During a series of exposures, differential tracking causes back- 
ground objects to shift away from the field of view limited by the 
22" wide strip mask used in polarimetric mode. However, since 
individual exposures were short (between 30 s and 100 s), star 
trailing was limited always to less than 1 pixel size, compared to 
a typical seeing of 4 pixels. 

The results of our test are illustrated in Fig. [4] which shows 
that for the background star N' c {a) tends to zero for a <: 
20 pixels, while for 133P/Elst-Pizarro N' c {a) converges to a posi- 
tive constant value. In Fig. [4] the constant used to disentangle the 
An(o) and Ac(fl) profiles was calculated by interpolating with a 
constant term the N' iol {a) profiles (shown in the middle panel) 
between 22 and 28 pixels. For both the background star and 
133P/Elst-Pizarro, all profiles were normalised imposing that the 
fluxes integrated within a 15 pixel (=3.75") aperture are equal to 
1. 

The most critical issue is that the coma region (if a coma 
is present at all) is contaminated by the tail contribution, which 
so far we have implicitly neglected. Therefore, we transformed 
our images into polar coordinates (a, 9) and removed the regions 
within those azimuth ranges that were clearly contaminated by 
tail(s) or background sources. The flux pertaining to a certain 
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annulus [a, a + da] was then calculated by integrating the pixel 
values at the various 9 values contaminated by netiher tail nor 
background sources, multiplied by a factor > 1 to account for 
the image trimming. The errorbars in the Nl, profiles were es- 
timated by associating with the measured flux the standard de- 
viation of the distribution of the fluxes at the various 9 values. 
Finally, we applied an algorithm similar to that described in the 
case of cartesian coordinates. The coma contribution was mea- 
sured in circular regions between 4500 and 9000 km from the 
nucleus photo-centre. 

The results of this analysis, shown in Fig. [5] are consis- 
tent with there being a coma with Afp £ 1 cm detected at a 
2-5 <t level. This result could well be ascribed to stray light in 
the instrument (which may mimic a diffuse halo). We there- 
fore repeated the same analysis on the images of a number of 
background stars at various observing epochs. We found that 
the ratios Nc(a) / N t ot(a) were generally substantially higher for 
133P/Elst-Pizarro than for background objects. In conclusion, 
our data exclude the presence of a coma with Afp k 3 cm, but 
certainly do not allow us to rule out the possibility that a faint 
coma, with Afp of the order of 1 cm or smaller, exists. Data 
shown in Fig. [5] marginally suggest a change of the colour in- 
dex, but reaching a firm conclusion in that respect would require 
higher signal-to-noise ratio data. 

This v alue compares to a dust p roduction rate of the order of 
100 gs~' (Boehnh ardt et alj |2008). which would be one of the 
lowest level ever measured for a comet, although lasting most 
likely for several months around perihelion. Despite its large un- 
certainty (order of one magnitude ), this du s t pro duction rate is 
still higher than that obtained bv lHsieh et alJ d2004l) . which could 
be caused by temporal variability in the nucleus activity and/or 
measurement and modeling errors (in both datasets). The most 
logical explanation of the almost absent coma and the narrow 
dust tails is low nucleus activity and the small terminal expan- 
sion velocity of the dust grains after release by the nucleus of 
133P/Elst-Pizarro. Small dust expansion velocities can be con- 
cluded from estimations of the tails. Hence, the dust coma is 
weak and confined very much to the near-nucleus region, which 
makes it difficult to detect in the atmospheric seeing disk of the 
latter. 

3.3. Nucleus photometry 

For each observing night, our data set typically consists of an 
acquisition image in R, a series of 12 to 32 polarimetric images 
in V, followed by a series of 12 to 32 polarimetric images in 
R, and a series of two photometric images in R and two in V. 
Images were taken approximately two or three minutes apart, 
except for the photometric series, which was occasionally taken 
ten to thirty minutes before or after the polarimetric series. For 
each epoch, the set of photometric data consists of the acquisi- 
tion image of the polarimetric series in R, the averaged frames 
obtained from all polarimetric images in R and V, and images 
obtained in R and V with no polarimetric optics. 

For all images, we used aperture photometry, adopting an 
aperture radius of 12 pixels (= 3")- Our nucleus brightness mea- 
surements are contaminated by the emission of a tail and (pos- 
sibly) a coma. Tail emission was roughly taken into account by 
subtracting, from the flux integrated in a circular aperture of ra- 
dius a (in pixels), the amount corresponding to Af = 2.5 x 10 -8 
per pixel, multiplied by an area of 4 (a - 2) pixels, both for the 
R and the V filter, for all observing epochs apart from the one 
before perihelion, when the tail was very faint. For the adopted 
12 pixel aperture, this corresponds typically to a 0.1 mag correc- 
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Fig. 6. Photometry in the R (red circles) and V (blue squares) 
filters for 133P/Elst-Pizarro as function of time, obtained on 22 
May 2007. Filled symbols refer to the photometry obtained from 
the polarimetric images, empty symbols refer to the photom- 
etry obtained from the images with no polarimetric optics in. 
Asterisks show the photometry obtained from the averaged po- 
larimetric images. Point at t — is obtained from the acquisition 
image of the polarimetric series. 

tion. Coma contribution was subtracted using the results of the 
previous section. For a 12 pixel aperture, this corresponds to a 
correction generally of the order of 0.1 mag. Finally, the appar- 
ent magnitudes were converted into normalised magnitude using 



H m (a) = m-5 log(rA) 



(8) 



where r and A (the heliocentric and the geocentric distances), are 
expressed in AU, and a is the phase-angle. 

For each observing epoch, we finally obtained a light curve 
with a typical time baseline of the order of 1 h 15 m to 1 h 30 m. 
At all the observing epochs, we could clearly note that the ob- 
ject photometry varyed within a short timescale (< 1 h). The 
most extreme case, shown in Fig. |6l is night 22 May 2007, 
during which we measured an amplitude of about 0.5 mag and 
0.6 mag in R and V, respectively, to be compared to a variabil- 
ity of S; 0.02 mag measured for a nearby background star. In the 
observations obtained on 19 September 2007, which have the 
shortest time baseline of our dataset (1.1 h), 133P/Elst-Pizarro 
exhibited the smallest brightness variability $ 0.2 mag in both 
filters. While in some cases, the variability of 133P/Elst-Pizarro 
may be explained by the close presence of background objects 
or by changes in the sky transparency, in m ost of the cases the 
observed variability is intrinsic to the object. iHsieh et al] {2004) 
estimated a light curve amplitude of the order of 0.4 magnitude 
in R filter, with a 3.471 h period. Earlier observation s reported 
substantially larger variability (see IHsieh et al1l2004L and refer- 
ences therein), which may indicate changes of the nucleus aspect 
angle along the orbit (although some results reported previously 
may be affected by poor calibration). 

We note that our observational dataset may possibly be 
suited to a detailed period analysis based on differential pho- 
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Table 2. Photometry of comet 133P/Elst-Pizarro. 



Date 


Phase 








H R {a) (mag) 




Hy (a 


) (mag) 


(yyyy mm dd) 


o 


(AU) 


(AU) 


acq} ' 


aver, pol. 1 ' 


photom. p; 


aver, pol. 1 ' 


photom. p ' 


2007 05 22 


19.69 


2.6467 


2.0129 


16.69 


16.76 


16.53 


17.26 


16.82 


2007 07 13 


3.33 


2.6424 


1.6329 




15.76 


15.62 


16.11 


15.94 


2007 07 17 


1.66 


2.6429 


1.6282 


15.57 


15.75 


15.53 


16.04 


15.88 


2007 08 01 


4.95 


2.6455 


1.6466 


15.68 


15.86 


15.65 


16.22 


16.13 


2007 08 05 


6.68 


2.6465 


1.6619 


15.91 


15.69 


15.92 


16.26 


16.40 


2007 08 09 


8.20 


2.6475 


1.2794 


15.95 


15.82 


15.91 


16.27 


16.35 


2007 08 30 


15.32 


2.6548 


1.8342 


(16.19) 


(16.34) 


(16.51) 


(16.33) 


(16.82) 


2007 09 19 


19.61 


2.6646 


2.0521 


16.46 


16.25 


16.36 


17.00 


16.94 


2007 09 24 


20.29 


2.6674 


2.1130 


16.42 




16.74 


16.63 


17.22 




0.00 








15.50 ±0.05 




15.88 


±0.07 



(1) Helio-centric distance to the object. (2) Geo-centric distance to the object. (3) Magnitude measured in the acquisition images. (4) Magnitude 
measured in the averaged polarimetric frames. (5) Magnitude measured in the frames obtained with no polarimetric optics in. The last row 
refers to the brighntess at phase-angle 0° extrapolated with a linear best-fit. Blanks refer to frames where the image of 133P/Elst-Pizarro was 
contaminated by background objects. Values given in parenthesis should be considered as lower limits because obtained during a cloudy night. 
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Fig. 7. Photometry in the R (red circles) and V (blue squares) 
filters for 133P/Elst-Pizarro as function of the phase angle. Solid 
lines show the corresponding linear best-fits. Points that were not 
used to calculate the best-fit are represented with empty symbols. 



tometry, which is beyond the scope of this paper. The aim of 
our photometric measurements is to determine, by extrapolating 
of the brightness measured at various phase angles, an estimate 
of the absolute brightness of the object nucleus at phase angle 
0°. With this in mind, the uncertainty introduced by using sim- 
ple aperture photometry (as opposed to the PSF fitting method), 
and by a less than optimal photometric characterization of the 
observing nights, is negligible compared to the variability of the 
object due to its rotation. The major contribution to the errorbars 
comes from the uncertainty in the zeropoints, and in both the tail 
and coma contributions, the sum of which was finally estimated 
to be 0. 1 mag for all points. 

Our results are given in Tableland plotted in Fig. [7] The 
slopes of the brightness curves are 0.052 + 0.004 mag/deg in R, 
and 0.055 + 0.005 mag/deg in V, indicating a phase darkening 



coefficient in the range typical of cometary nuclei. The extrapo- 
lated average brightnesses at phase angle = 0° are 15.50 + 0.05 
and 15.88 + 0.07 in the R and V filters, respectively, resulting 
in an average colour index V — R = 0.38, which is equivalent 
to an average spectral gradient of +2%/100nm in the V and R 
wavelength range, i.e., the intrinsic color of 133P/Elst-Pizarro 
appears to be neutral, and a flat solar reflectance spectrum in the 
visible is expected, at least beyond about 500 nm. 

The absolute brightness of 133P/Elst-Pizarro in R and V fil- 
ter clearly increases towards zero phase angle, but no indication 
of an opposition brightening is found. Deviations from linear- 
ity may be caused by measureme nt errors and (mos tly) rotation 
variations of the elongated body (Hsie h et al.ll2004 give a min- 
imum axis ratio of 1.45 and a rotation period of 3.471 h). The 
absolute brightnesses at zero phase-angle correspond to an aver- 
age equivalent radius of 1.6 km (range is 1.64-1. 72 km in R to 
1.47-1. 57 km in V), assuming the 0.07 and 0.06 albedo values 
(in R and V, respectively) obtained from our polarimetric mea- 
surements. 

We finally compare our results with those of previous stud- 
ies, noting that previous pho tometric measurements were ob- 
tained more than fo ur months (iBoehnhardt et al.ll 19971) and more 
than seven months dHsieh et al.l l2004) after perihelion in 1996 
and 2002, respectively, while the measurements presented in this 
work cover the orbit arc from about one month before to three 
months after perihelion in 2007. 

iHsieh et all (120041) and Uewitt et all (120071) found a phase 
darkening coefficient /Jd = 0.044 +0.007 mag/deg and an average 
V-R color of 0.42+0.03, which are fully consistent with our own 
measurements. We note that colour indices obtained from quasi- 
simultaneous images in R and V vary bet ween 0.29 and 0.48 , 
which is marginally larger than reported by iHsieh et al.l (12004). 
who measured V-R varying between 0.35 and 0.49. 

The nucleus size has been previously estimated (neglect- 
ing possible contamination from coma and tails) to be about 
2.5 ± 0.1 km assuming an albedo of 0.04, which is also consis- 
tent with our value s when applying identical albedo parameters. 
Hsieh et al. 1 (120091) give an equivalent radius of the nucleus of 
133P/Elst-Pizarro of 1.9 ±0.3 km for a geometric albedo of 0.05. 

We note that the variation amplitudes of 0.5 to 0.6 mag seen 
in our photometric data over time intervals between 1 h and 1 .5 h, 
compare to a minimum aspect ratio (ratio of the long to short 
axes lengths for a prolate ellipsoid rotating about its small axis) 
of about 1.6 to 1.7. These values are slightly larger than that 
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Fig. 8. The measured Pq values as a function of the phase angle 
in the R (red circles) and V (blue squares) filters. The solid lines 
show the best-fit relation obtained with a second order polyno- 
mial expansion constrained to pass through (0,0). 
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9. Distribution of the Py values expressed in units of error- 
The vertical dotted lines refer to the 3-cr detection level. 



measured by (Hsieh et al. 2004), which could be due to the dif- 
ferent aspect angles of the varying nucleus cross-section along 
the orbit of the object. 

3.4. Imaging polarimetry 

Polarimetry was calc ulat ed using the method outlined, e.g., by 
Bagn ulo et al.l( 2006) and Bagnulo et al.l d2008l) . Null parameters 
(see iBagnulo et al.l l2009}) were also systematically calculated to 
check the reliability of the results. The final reduced Stokes pa- 
rameters Pq and Pjj were obtained by adopting as a reference di- 
rection the perpendicular to the gr eat circle passin g through the 
object and the Sun, using Eq. (5) of Ba gnulo et al.l d2006l) . In this 
way, Pq represents the flux perpendicular to the plane defined by 
the Sun, Object, and Earth (the scattering plane) minus the flux 
parallel to that plane, divided by the sum of these fluxes. Our fi- 
nal polarimetric measurements are reported in Table [3] Figure [8] 
shows Pq results as a function of the phase angle in both filters. 

For symmetry reasons, Py values are always expected to be 
zero, and inspection of their values allows us to perform an indi- 
rect quality check of the Pq values. Figure|9]shows the distribu- 
tion of the Py values in the R and V filters expressed in errorbar 
units. Since the distribution appears to be peaked at about zero, 
and none of the Py values exceed a 3-cr detection, this test repre- 
sents a positive quality check of our polarimetric measurements 
of 133P/Elst-Pizarro. 

3.4.1. Polarimetry of the tail 

We attempted to measure the polarization of the tail, but, not 
unexpectedly, the low signal-to-noise ratio prevented us from 
obtaining accurate measurements. The highest precision was 
reached in the images obtained on 13 July 2007 and 17 July 
2007, where Pq and Py could be measured with a ~ 3 % er- 
rorbar. In the R filter, we integrated the signals over an area 
of 110 and 130 pixels, respectively, and we obtained Pq - 



-1.60 ± 3.0% and Py = -0.5 ± 3.0% on 13 July 2007, and 
P Q = -1.1 + 2.5% and Py = 0.4 + 2.5% on 17 July 2007. 
On the image obtained on 5 August 2007, we integrated over a 
40 pixel area of the tail, and we obtained Pq = -1.0 + 3.8 % 
Py - 1.6 + 3.8 %. In all cases, the null parameters were consis- 
tent with zero at the 1.5 cr level. 

3.4.2. Identifying the polarization of the nucleus 

We now discuss whether our measurements can be considered 
representative of the nucleus of 133P/Elst-Pizarro, and to what 
extent they are contaminated by the coma and tail of the object. 

We used the results obtained in the polar coordinate system 
to estimate the true fluxed produced by the nucleus of 133P/Elst- 
Pizarro within the aperture used for our polarimetric measure- 
ments (see col. 8 of TableO, and we found that for the apertures 
selected for the polarimetric measurements, the contributions 
from the coma were between 0.05 and 0.1 mag, while the contri- 
bution from tail was generally 0.02-0.03 mag. We can therefore 
assume that, for the aperture values used for polarimetry, at least 
about 90 % of the measured flux is related to the nucleus con- 
tribution. Assuming that the flux contribution due to coma and 
tail is entirely produced by a dust scattering mechanism rather 
than gas emission, we conclude that the radiation due to coma 
and tails is either parallel to the light reflected by the nucleus, or 
perpendicular to it. This hypothesis is supported by all of our Py 
measurements being consistent with zero. Following Eq. (1) of 
Bagnu lo et al . (2008), we can thus write the nucleus polarization 



P 1 * = UP™ S -0.1P£ +T , (9) 

where fg bs is the total observed polarization, and P^ +T is the 
fraction of linear polarization of the light scattered by coma and 
tail. For instance, if coma and tail are not polarized (Pq +t =0), 
the nucleus polarization is underestimated (in absolute value) by 
about 10 %; if coma and tails are polarized at a 3 % level parallel 
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Table 3. Polarization measurements of 133P/Elst-Pizarro. 



Date 


TJT<1) 


Exp 


PA <2) 


Phase Angle 


Ci) 

t — t 


Filter 


Ap. <4) 


Pq 


Pu 


mm r\r\ 

yyyy 


n n ■ mm 
1111.111111 


W 


*. ) 


K ) 








\ 10 ) 




2007 05 22 


09:44 


960 


256.160 


19.690 


-37.66 


R 


7.0 


0.96 ±0.59 


0.10 ±0.59 


2007 05 22 


09:17 


960 


256.162 


19.693 


-37.68 


V 


7.0 


0.85 ±0.55 


-0.01 ±0.56 


2007 07 13 


07:43 


960 


247.877 


3.331 


+ 14.25 


R 


5.0 


-0.75 ±0.16 


-0.39 ±0.16 


2007 07 13 


07:13 


960 


247.905 


3.340 


+ 14.23 


V 


5.0 


-1.27 ±0.15 


-0.10 ±0.16 


2007 07 17 


07:05 


960 


237.251 


1.646 


+ 18.22 


R 


7.0 


-0.53 ±0.16 


-0.16±0.16 


2007 07 17 


06:37 


960 


237.360 


1.654 


+ 18.21 


V 


5.5 


-0.53 ±0.15 


-0.04 ±0.1 5 


2007 08 01 


02:47 


960 


85.710 


4.969 


+33.05 


R 


5.5 


-1.17 ±0.31 


0.02 ±0.29 


2007 08 01 


02:19 


960 


85.720 


4.960 


+33.03 


V 


7.0 


-0.50 ±0.39 


0.09 ±0.39 


2007 08 05 


04:58 


1280 


84.169 


6.671 


+37.14 


R 


5.0 


-1.28 ±0.14 


0.26 ±0.15 


2007 08 05 


04:27 


1280 


84.174 


6.662 


+37.12 


V 


6.0 


-1.48 ±0.15 


-0.14 ±0.15 


2007 08 09 


01:09 


960 


83.378 


8.215 


+40.98 


R 


5.0 


-1.22 ±0.17 


0.12 ±0.17 


2007 08 09 


00:42 


960 


83.381 


8.207 


+40.96 


V 


5.0 


-1.53 ±0.18 


0.02 ±0.18 


2007 08 30 


00:43 


1280 


81.907 


15.330 


+61.96 


R 


6.0 


-0.28 ±0.46 


0.03 ±0.48 


2007 08 30 


00:13 


1280 


81.908 


15.324 


+61.94 


V 


6.0 


-0.51 ±0.29 


-0.02 ±0.30 


2007 09 19 


02:16 


1200 


80.962 


19.611 


+82.02 


R 


6.0 


0.33 ±0.37 


-0.32 ±0.36 


2007 09 19 


01:48 


1200 


80.963 


19.608 


+82.00 


V 


5.0 


0.99 ±0.39 


0.08 ±0.37 


2007 09 24 


01:01 


1280 


80.655 


20.288 


+86.97 


R 


7.0 


1.63 ±1.13 


0.27 ±1.17 


2007 09 24 


00:28 


1280 


80.656 


20.285 


+86.95 


V 


6.5 


1.94 ±0.90 


-0.64 ±0.88 



(1) Date and UT refers to the middle of the exposure series. (2) Position angle of the scattering plane. (3) Time from perihelion. (4) Aperture 
used for flux extraction, expressed in pixels (1 pixel = 0.25"). 

to the scattering plane (Pq +t = -3 %), the observed polarization 
minimum of about — 1.5 % would correspond, for the nucleus, 
to ~ -1.35%; if coma and tails are polarized at a 3 % level 
perpendicular to the scattering plane (for the case of Fresnel- 
reflection mechanism), the observed polarization minimum will 
be Pq = -1.95 %. We note that our direct measurements pre- 
sented in Sect. 13.4.11 although affected by large errorbars, seem 
more consistent with a negative value of the tail polarization than 
a positive one, and infer that the fraction of tail linear polariza- 
tion is very similar to the fraction of polarization from the nu- 
cleus plus tail. Our dataset does not allow us to go beyond these 
qualitative statements, yet our conclusion is that the observed 
polarization is representative of the polarization of the nucleus 
of 133P/Elst-Pizarro. The following analysis will be based on 
this conclusion. 



3.4.3. Characterization of the polarimetric curve 

As for the large number of atmosphere-less solar system bodies, 
the nucleus of the 133P/Elst-Pizarro exhibits the phenomenon 
of negative polarization: at small phase angles, the electric field 
vector component parallel to the scattering plane predominates 
over the perpendicular component, in contrast to what is ex- 
pected from the simple single Rayleigh-scattering or Fresnel- 
reflection model. This phenomenon is generally explain ed in 
terms of coherent backscattering (e.g.. iMuinonenl 120041) . For 
small bodies of the solar system, the polarization reaches a min- 
imum value generally between phase angle 7°and 10°, and be- 
comes positive at p hase angles of between 17° and 22° (e.g., 
Penttila et al. 2005). In conclusion, for those small bodies for 
which it is possible to measure the behaviour of the polarization 
for an extended range of phase angles, three important charac- 
teristics can be measured well, and used to perform an empirical 
classification: the minimum polarization and its corresponding 



phase angle, the slope of the polarimetric curve in the linear 
part beyond polarization minimum, and the inversion angle at 
which the polarization changes from being parallel to the scat- 
tering plane and becomes, at larger phase angles, perpendicular 
to the scattering plane. Following these guidelines, we perform 
an analysis of our polarimetric data for the nucleus of 133/Elst- 
Pizarro. 

Inspection of Fig. [8] shows that measurements in both the R 
and V filters agree within the errors (except for a point at phase 
angle 3°). Nevertheless, we perform our basic analysis based on 
the measurements obtained in the R and V filters separately. 

The best-fit relation obtained using a second order polyno- 
mial (constrained to pass through the origin) shows that in the R 
filter, the minimum is reached at phase angle a™ 11 = 8.8° + 0.6°, 
for a minimum polarization value of P™ n (R) = -1 .30 + 0.20 %. 
In the V filter, the minimum is reached at phase angle = 
8.5° + 0.5°, with a minimum polarization value of P™ n (V) = 
— 1.57 ± 0.20% . The positive crossover, when the observed po- 
larization changes sign from negative to positive, is at phase an- 
gle a R (P Q = 0) = 17.6° + 2.1° in the R filter, and at phase angle 
a v (P Q = 0) = 17.0° + 1 .6° in the V filter. 

Using a second order polynomial, the slope of the polariza- 
tion depends clearly on the phase angle in which it is calculated. 
With reference to the the crossover point, the polarization slopes 
are 

l dPg(R) \ = 0.30 + 0.02 %/deg 
\ da lp a =o 

and 

( f | = 0.37 + 0.02 %/deg , 
\ da )p =o 



S. Bagnulo et al.: Polarimetry and photometry of 133P/Elst-Pizarro 



11 



for the R and the V filter, respectively. We note however that 
these estimates are obtained based on the assumption that the 
measurements may be fitted with a second order polynomial. 

3.4.4. Relationships between polarization and albedo 

For asteroids, two empirical relationships have been found be- 
tween the polarization characteristics and the albedo, namely, 
both slope and polarization minimum seem to be related to the 
geometric albed o of the body by means of a si mple expression 
as given, e.g., bv lLupishko & Mohamedld 19961) 

log(p) = -0.98 log ^) -1.73 (10) 
and 

log(/>) = -1.22 log (p™")- 0.92, (11) 

where p is the albedo, dPg/da is the derivative of the observed 
linear polarization with respect to phase angle, and P™ n is the 
minimum value of the polarization. Using Eq. ([Tol l, from our 
data we obtain for the albedo the values pr - 0.061 + 0.004 in 
the R filter, and p v = 0.049+0.003 in the V filter. Using Eq. {TT), 
we obtain p R = 0.087 + 0.017 and p v = 0.069 + 0.011. Their 
average values correspond to about pr = 0.07 and pv = 0.06, 
uncertainties being about 0.01. 

The average albedo value obtained from our polarization 
measurements agrees remarkably well with that determined via 
the classical method of combined visible and thermal flux mea- 
surements of 133P/Elst-Piza rro, i.e., /?« = 0.05 + 0.02 as pub- 
lished by Hsie h et al.l (|2009). However, it remains to be shown 
that the albedo-polarization relationships are indeed applicable 
to our target since they were verified and calibrated using aster- 
oid data for different taxonomic types than those that we have 
found for 133P/Elst-Pizarro (see Sect.©. 

4. Discussion 

There are two main aspects to the properties of 133P/Elst- 
Pizarro. On the one hand, it is located in the asteroid belt with 
an orbit similar to the Themis collision family or a subgroup of 
it, the Beagle family, and over a large part of its orbital revolu- 
tion it appears asteroid-like. On the other hand, around and after 
perihelion it displays activity producing a dust tail, as comets 
do. In the following, we discuss whether our new polarimetric 
measurements could help us to identify whether it is an asteroid 
(with cometary activity) or a comet (scattered into the asteroid 
belt). 

We first compare the polarimetric curve of the nucleus of 
1 3 3P/Elst-Pizarro with th at of the nucleus of comet 2P/Encke 
dBoehnhardt et al.l [2008) by making use of the scattering 
parametrisation provided in iMuinonen et all d2009allbl) . By us- 
ing lines of electric dipoles with an inverse-gamma distri- 
bution for the line lengths, and by estimating the parame- 
ters with a Markov-Chain M onte-Carlo (MCMC) method (see, 
e.g.. IMuinonen et aLll2009bl) . we have computed, for 133P/Elst- 
Pizarro and 2P/Encke, physical interference models. The eight 
parameters of the full model are: the exponents or shape param- 
eters S\ and S2 and the scale parameters y\ and y2 of the inverse- 
gamma distributions for the two polarization states (subscripts 1 
and 2 are for positive and negative states, respectively); the com- 
plex amplitude Z and interdipole distance kd^ (k - 2n/A is the 
wave number and dj, the physical distance) for the longitudinal 



Table 4. MCMC parameters for 133P/Elst-Pizarro, 2P/Encke, 
and (1) Ceres. S X(Z) and 3(Z) are the real and imaginary parts 
of the complex amplitude Z, respectively. 
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electric-dipole contribution; and the normalized weight w of the 
Rayleigh-like contribution to the polarization curve. 

For 133P/Elst-Pizarro and 2P/Encke, reduced six parameter 
models already suffice to explain the polarimetric observation. 
We assume that s\ — S2 — s and fix kdi - n. Least squares 
fitting analysis and MCMC sampling provide the polarization 
curve corresponding to the best-fit relations (solid lines), and in 
the 3 cr error envelopes (dashed lines) depicted in Fig. [10] Best- 
fit model parameters and their 3 <x errors are given in Table [4] 
Among the parameters, the weights w, exponents s and scale 
parameters yi and y^ are of importance for the present study. 
The weights w obtained realistic values w > 0.5, which, in the 
physical model, indicate a positive contribution from point-like 
Rayleigh scatterers in addition to the contributions from the lines 
of dipoles. The exponents s and the scale parameters y\ and yi 
are loosely constrained by the polarimetric data in the current 
case of considerable observational errors. 

In addition to 2P/Encke and 133P/Elst-Pizarro, the physical 
model has so far b een applied only to the polarimetric data of 
asteroid (1) Ceres (IMuinonen et al.ll2009bl) . The corresponding 
best-fit model parameters are given in Table [4] 

We conclude that it is too early to reach definitive conclu- 
sions using MCMC models, beyond that the polarization char- 
acteristics of 133P/Elst-Pizarro are closer to those of (1) Ceres 
than those of 2P/Encke. The polarimetric slope at the inversion 
angle is larger for 133P/Elst-Pizarro than for (1) Ceres. All three 
weights w can overlap within their 3-cr error domains and the 
same is true for the scale parameter yi. Mainly in terms of the 
scale parameter yi, 2P/Encke differs from the others. 

A more statistically meaningful approach is offered by a 
straight comparison of the photometric and polarimetric prop- 
erties of 133P/Elst-Pizarro with those of different types of aster- 
oi ds, cometary dust, Cen taurs and Kuiper-Belt objects. Table 2 
of Bo ehnhardt et al.l d2008l) lists the typical numerical values for 
albedo, slope of the photometric phase function, spectral gra- 
dient, minimum polarization and the related phase angle, slope 
of the polarization phase function, inversion angle, and spec- 
tral gradient of the polarization for different classes of small 
bodies of the solar system. A comparison with the values of 
133P/Elst-Pizarro derived in this work allows us to conclude 
that 133P/Elst-Pizarro does not exhibit the typical properties of 
Kuiper-Belt objects and Centaurs, C, S, E and M type asteroids, 
cometary dust, while there is a rather good agreement in all pa- 
rameters with those of F-type asteroids. 

F-type asteroids are considered to be primitive dBus et al.l 

2002), and t hey are sometimes cla imed to be related to comet- 
like objects dWeissman et al.f 2002). A weakness in this compar- 
ison is the meager knowledge of polarimetric parameters of cer- 
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Fig. 10. Polarimetric observations of the nucleus of comet 2P/Encke (left panel) and of 133P/Elst-Pizarro (right panel), modeled 
using MCMC for the empirical linear-exponential model. Blue squares refer to V Bessell filter, and red circles to R Bessell filter. 
For 2P/Encke, red and blue asterisks refer to narrow band filters centred about 834 and 485 nm, respectively. The best-fit models are 
shown with solid line, while the 3 <x error envelopes are shown with dashed lines. 



tain types of small bodies, such as cometary nuclei, Centaurs, 
and also F-type asteroids themselves (although the polarimetric 
dataset for F-type asteroids is substantially larger than that for 
Centaurs and cometary nuclei). However, if the parameter sim- 
ilarity between 133P/Elst-Pizarro and F-type asteroids is to be 
seen as an indicator of the similarity of these objects, one can 
predict that more objects with cometary activity will be found 
among asteroids with F-type taxonomy. 



Licandro et al. | (120071 assigned 133P/Elst-Pizarro a C- or 
B-type classification taxonomy classes to which F-type as- 
teroids are closely related to (the water absorption feature 
around 3/mi is missing and differences in the continuum be- 
low 0.4 jum exist for F-type objects). Other notable, pecu- 
liar C-class objects include 10 7P/Wilson-Harrington (C- or F- 
tvpe:lTholen&Barucci|ll989h. 3200 Phaefhon (B- or F-type; 
Thol en&Baruccilll989l:lLicandro etalJl2007bt). and 200 5 UD 
(B- or F-type: Uewitt&Hsieh Il2006t iKinoshita et alj|2007l) . The 
former two objects at least have displayed cometary activity 
in the past. An F-type 133P/Elst-Pizarro would thus have good 
company with respect to cometary activity. 

Comet migration into the asteroid be lt does not seem t o 
be a very efficient process dynamically dLevison et al. 1 12006). 
and the "snowline", i.e., the distance to the Sun at which wa- 
ter vapour in the protoplanetary disk condenses and becomes 
accreted in forming planetesimal, falls into the (outer) aster- 
oid belt. Therefore, one may be inclined to accept 133P/Elst- 
Pizarro as a "child of the inner solar system", and not of the cold 
outskirts where cometary nuclei originate. Our observational re- 
sults are certainly compatible with this scenario, although it im- 
plies that cometary activity can also occur in asteroids. The dust 
activity in 133P/Elst-Pizarro and the other so-called 'main-belt 
comets' require a driving mechanism. Ice sublimation is consid- 
ered a possibility, but is not as yet supported by observational 
evidence. 



5. Conclusions 

With the FORS1 instrument of the 8 m ESO VLT, we have car- 
ried out imaging observations in polarimetric mode of the object 
133P/Elst-Pizarro at nine observing epochs. The first observa- 
tion was performed in May 2007, about one month before peri- 
helion, and the remaining eight from July to September 2007, up 
to about three months after perihelion. These observations cover 
the phase-angle range of 0°-20°. 

Our images have detected one or two tails close to perihe- 
lion, and, to a marginal level of evidence, the presence of a coma. 
We have then performed an analysis of the dust production. The 
dust release of 133P/Elst-Pizarro may have started about one 
month before perihelion passage and continued at least until 
about three months after perihelion. The onset of dust activity 
in 2007 happened earlier than concluded for the 1 996 appari- 
tion o f the comet bv lPravec & Sekanina 1 d 19961) and Boehnhardt 
(1996), although, admittedly, the latter analyses could not re- 
ally constrain dust release before perihelion passage. However, 
the dust release in 133P/Elst-Pizarro is clearly repetitive and it 
starts and extends around perihelion passage. The dust grains 
that dominate the optical appearance of the tail were of microm- 
eter size and larger, and were released at a relatively low speed. 
The appearance of the tail in 2007 resembles that seen in 1996, 
although the two tail patterns refer to different emission peri- 
ods with respect to perihelion passage. The 1996 tail appearance 
supports the scenario of dust production that continued for up to 
half a year after p erihelion passage. Considering the results of 
iHsieh et all d2004l) . the active period would even extend to about 
a year after perihelion. The dust tails seen in 133P/Elst-Pizarro 
during the 1995, 2002, and 2007 apparitions support the con- 
cept that the activity of 133P/Elst-Pizarro is recurrent and that 
it extends over the same arc of the orbit (from about 1 month 
before perihelion to at least 1 year after perihelion) . A possible 
yet unproven scenario, described in greater detail by Hsieh et al. 
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( 20041) and implicit l y assu med also in earlier publicat ions by 
iPravec & Sekaninal d 19961) and lBoehnhardt et all (Il997l) . is that 
a main region of local activity is switched on - possibly due to 
solar illumination - shortly before perihelion and continues to be 
active along the orbit over a year or so. The repetitive behaviour 
of the activity would suggest a rather stable orientation of the ro- 
tation axis over the past three apparitions of 133P/Elst-Pizarro. 
The activity level of the nucleus is not very high and produces 
micron-size dust at low speed. The true release mechanism could 
not be constrained from our observations. 

We have used our observations to characterize the polari- 
metric and photometric behaviour of the object as a function 
of its phase-angle, and we have compared it to that of the nu- 
cleus of comet 2P/Encke, and other small bodies of the solar- 
system. In particular, the comparison with the polarimetric curve 
of 2P/Encke was performed by adopting a physical classification 
tool for the polarization phase curves, but our results are not con- 
clusive because current lack of a statistically meaningful sample 
of objects that have been analysed in a similar fashion. However, 
a direct comparison of the polarimetric and photometric curves 
identified several similarities with the observational parameters 
of F-type asteroids. 
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